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Abstract. We study the quantum mechanical motion in the x 2 y 2 potentials with 
n = 2,3, which arise in the spatially homogeneous limit of the Yang-Mills (YM) 
equations. These systems show strong stochasticity in the classical limit (h = 0) 
and exhibit a quantum mechanical confinement feature. We calculate the partition 
function Z(t) going beyond the Thomas-Fermi (TF) approximation by means of the 
semiclassical expansion using the Wigner-Kirkwood (WK) method. We derive a novel 
compact form of the differential equation for the WK function. After separating the 
motion in the channels of the equipotential surface from the motion in the central 
region, we show that the leading higher-order corrections to the TF term vanish up to 
eighth order in h, if we treat the quantum motion in the hyperbolic channels correctly 
by adiabatic separation of the degrees of freedom. Finally, we obtain an asymptotic 
expansion of the partition function in terms of the parameter g 2 h 4 t 3 . 
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1. Introduction 

The discovery of the chaoticity of the classical Yang-Mills (YM) equations pQ (see [2] 
for a review) has attracted broad interest to the system of two (three) coupled quartic 
oscillators with the potential x 2 y 2 (x 2 y 2 + y 2 z 2 + z 2 x 2 ), where x,y (z) are functions of 
time t. These systems are the simplest limiting cases for the homogeneous YM equations 
(sometimes called YM classical mechanics) with n = 2 (n = 3) degrees of freedom, 
respectively. The x 2 y 2 model, the central object of the present paper, exhibits a rich 
chaotic behavior despite its extreme simplicity. Not surprisingly, this model has been 
encountered in various fields of science, including chemistry, astronomy, astrophysics, 
and cosmology (chaotic inflation). 

Quantum mechanically, this model (YM quantum mechanics, YMQM), despite 
possessing a logarithmically divergent volume of energetically accessible phase space 
[3 E] (we set m — 1 throughout) $ 



can be shown to have a discrete spectrum Ej- Physically, it is clear why this is 
so: Quantum fluctuations, e.g. zero-point fluctuations, forbid the "particle" to escape 
along the x or y axis where the potential energy vanishes. The system is thus confined 
to a finite volume, and this implies the discreteness of the energy levels. Classically, 
of course, the particle can always escape along one of the axes without increasing its 
energy. 

In this article we calculate the partition function (heat kernel) Z(t) for the 
YMQM with the above given potentials beyond the well-known Thomas-Fermi (TF) 
approximation, which takes into account only the discreteness of the quantum 
mechanical phase space, but treats the Hamiltonian classically. There are several 
interesting articles IHj devoted to the approximate calculation of Z(t) in the TF 
approximation for these potentials. They are based on the adiabatic separation of 
the dependence of Z on x and y in the narrow channels of the equipotential surface 
xy = const.. The range of the integration over the coordinates x, y and the momenta 
PxiPy was divided into two regions: the central region (|x|, \y\ < Q) and the channels 
(Q < \%\,Q < \y\), which are governed by quite different physics: The system 
is essentially classical in the central region and intrinsically quantum mechanical in 
the channels performing oscillatory motion with rc-dependent frequency in the "fast" 
variable y (in the channels along the x-axis), but quasi-free motion in the "slow" variable 
x. 

The dependence on the artificial boundary Q dividing the central region from the 
channels disappears in the final answer for Z(t). Below we show that this property 
survives beyond the TF approximation up to the order h 8 . In the paper [Hj an alternative 
to this approach of calculating Z(t) was proposed, starting from the Yang-Mills-Higgs 

| This is in violation of Weil's famous theorem, which states that the average number N(E) of energy 
levels with energy less than E is asymptotically proportional to J Te>cIE' . 




(1) 
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quantum mechanics (YMHQM) and then passing to the limit v — 0, where v is the 
vacuum expectation value of the "Higgs field" defining the strength of the harmonic 
potential. 

Here we study the x 2 y 2 potential and postpone the investigation of YMHQM to 
a separate publication We here follow the method of separation of the domain of 
motion into two regions, the central square x, y G [-Q, Q] and the hyperbolic channels, 
introduced by Tomsovic [7j and Whelan jS] . The motion in the channels can be treated 
by adiabatic separation of the motion in the slow variable (in the direction of the channel) 
and the motion in the fast variable (perpendicular to the channel). We will go beyond 
the Thomas-Fermi approximation used in IE] by applying the well-known Wigner- 
Kirkwood (WK) method [TT) I12[ ITS] (see ^3] for a review of the WK approach). 

We already mentioned that the motion in the central region, where the variables 
x and y are treated on an equal footing, is quite different from the motion in the 
channels, where the motion in the perpendicular direction (for now taken as the variable 
y) performs quantum oscillations with x-dependent frequency, but the motion in the x 
variable is not treated as free as it was done in [7J |H]. We apply the WK method in 
the central region to both variables, but only to the slow variable x in the channels 
treating the motion in the fast y variable fully quantum mechanically. We improve the 
method of adiabatic separation introduced in |H] by taking into account the effect 
of the x-dependence of the oscillation frequency onto the motion in the x-direction. 
This is especially important, as we shall see, for the higher-order quantum corrections. 
We show that, up to eighth-order in % and in the leading terms in {tQ A )~ l <C 1, the 
calculated partition function Z(t) does not depend on the boundary Q dividing the two 
regions and the quantum corrections from both regions cancel. 

An interesting phenomenon occurs in the channel motion due to the presence of 
terms suppressed by powers of (tQ 4 ) -1 , which are independent of Q. Each term of 
this kind is negligible in comparison with the leading (but canceling) terms, but their 
series forms an asymptotic series in the variable g 2 h 4 t 3 , which does not involve Q. If we 
postulate that the entire Q-dependence of the partition function disappears, which we 
prove for the leading terms in (iQ 4 ) -1 , only this series will be left as the contribution 
from the channels. We may express this phenomenon as the "transmutation" of the 
small expansion parameter (tQ A )~ l governing the approach of adiabatic separation of 
variables into the small parameter g 2 h A t 3 characteristic of quantum mechanics. The 
quasiclassical expansion is shown to have the nature of an asymptotic series. 

In the next two sections we present the YMQM system and the WK method of 
calculating Z(t) beyond the TF approximation. 
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2. Yang-Mills-Higgs classical and quantum mechanics: The Thomas-Fermi 
approximation 

For spatially homogeneous fields (long wave length limit of the Yang-Mills field) the 
classical Hamiltonian for n = 2 is given by the expression 



l (J2 , „2X , 9 



2 



H = -(pI+pI) + »-xV. (2) 
The quantized counterpart of (J2J) is 

H = -%Vl y + • (3) 



A brief word on units: All quantities are given below in units of the energy E with 
dimensions [H] = 1, [t] = —1, [x], [y] = 1/4, [g] = 0, [h] = 3/4. The operator (jSJ) has a 
discrete spectrum jSJEl- The TF approximation to the heat kernel or partition function 
Z(t) = Tr[exp(— tH)} is the standard lowest-order semiclassical approximation valid for 
small M 3 ^ < 1. It is obtained by substituting the classical Hamiltonian for its quantum 
counterpart and replacing the trace of the heat kernel by the integral over the phase- 
space volume normalized by (2iTh)~ n , where 2n is the phase-space dimension. In other 
words, the TF approximation takes into account only the discreteness of the quantum 
mechanical phase space, but considers momenta and coordinates (in our case, the field 
amplitudes x and y) as commuting variables. This method was used in numerous papers 
(see e.g. jHISE])- For the calculation of the energy level density p(E) = dN(E)/dE at 
asymptotic energies, the TF approximation is a consistent approach since, as we shall see 
below, all corrections to the TF term are structures with factors h k t £ with k, £ positive 
integers. For the asymptotic energy level density p(E) or N(E) these corrections are 
negligible according the Karamata-Tauberian theorems E] relating the most singular 
part of Z(t) to the asymptotic level density: N(E) = J dEp(E) = L _1 (Z(t)/t) where 
denotes the inverse Laplace transform. 

For the Hamiltonian (jSJ) the naive TF approximation is divergent, because the 
classical phase space is infinite. As explaines in the Introduction, finite results for Z(t) 
and N(E) can be obtained by including certain quantum corrections to the channel 
motion by means of the method of adiabatic separation of the motion along and 
perpendicular to the channels [3 E]. We here give the expression for the partition 
function obtained in this improved TF approximation in our notation: 

1 ln-^- + 91n2 + C) , (4) 



u V^ghH^ 2 \ g 2 h A t 3 

where C is the Euler constant. We shall return to (@J) again below. Because we shall 
often encounter the pre-factor appearing in (j3J), we introduce the special symbol for it: 

K = (2ng 2 hH 3 )- 1 / 2 . (5) 
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3. Beyond the TF approximation: The Wigner-Kirkwood expansion 

For non asymptotic energies and also for the fluctuating part of the level density one 
needs to go beyond the TF approximation and calculate the quantum corrections to the 
TF term in Z(t). We remark that this problem is interesting not only from this practical 
point of view, but also because it provides insight into some problems of perturbation 
theory in quantum mechanics. Here we use the Wigner-Kirkwood (WK) method. Before 
we proceed, we emphasize the dual role of the variables x(t),y(t), z(t). They formally 
play the role of the coordinates, but at the same time they stand for the homogeneous 
gauge field amplitudes. The trace in the heat kernel may be taken with respect to any 
complete set of states. For the semiclassical expansion involving integration over the 
phase space for the quantum operators in the Wigner representation it is convenient to 
use the plane waves as a complete set: 

Z(t) = j—^ J n dxidpi e -^l n e' m e^/ h (6) 

with f = (xi,X2, • • • , x n ),p = (pi,p 2 , . . . ,p n )- The kinetic energy term in the exponent 
of exp(— tH) does not commute with the potential energy V(x); the WK expansion 
provides a convenient method of calculating the noncommuting terms. Following ^ 

we set 



tHJpr/h _ -tH, 



e e 



i ^l h W{f,p-t) = u(r,p-t) (7) 



where the function W(r, p;t) is to be determined. H(p,r) (as opposed to H) is the 
classical Hamiltonian (j2J). The function u(r,p;t) satisfies the Bloch equation: 

! + *« = <> (8, 

with the boundary condition 

\imu(f,p ] t)=e i ^ n , (9) 

corresponding to the initial condition W(r, p;0) = 1. From ((7|) and (jH)) we obtain an 
exact equation for W: 

2it 



dW h 



2 



dt 



A - t(AV) - ^(p • W) + t 2 {VVf 
a 



2 

+ T (ip- MVV) ■ V 

h 



W } (10) 

lb J 

where A = V 2 is the Laplacian. Next we expand W in powers of h: 



w = 52h k w k (n) 

k=0 

and equate the terms with the same power of h on both sides. We thus obtain a 
recurrence relation of differential equations for the Wk'- 
dW k 1 



k-2 



y [A - t(AV) + t 2 {VVf - 2tVV ■ V 
+ z£.[V-*(W)]W fc _i, (12) 
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with the initial conditions W k = for k < 0,Wq = 1. Since Z(t) and W are linearly 
related, the expansion of W in powers of h immediately translates into an expansion of 
Z(t) in powers of %. To obtain the term Zk(t), one needs to calculate W k from (fT2"J) and 
integrate over p and x: 



z k (t) 



Y[ dxidpiW k (r, p; t) exp 



-t 



P' 



+ V(x) 



(13) 



The expressions f!10|) and (|T2*j) may be written in more compact form, if we introduce 
the "covariant derivative" D = V — tVV^: 

<9iy ft 2 r„o -1 „. 1 r 

Y 



at 

or its recursive form 
dW k 



D 2 + -p-D 

n 



W 



(HD + ip 



V 



w 



(14) 



1 r 



D 2 W k „ 2 + 2ip-DW k . 



(15) 



dt 2 

We emphasize that the symbol p here denotes a classical phase-space variable and not 
an operator. As far as we know, the form (fPfj) of the equation (jlUj) has not previously 
been presented in the literature. 

In terms of the operator D, eq. (|14j) resembles a Fokker-Planck equation for 
W{r,p;t) with the diffusion constant a = h 2 and the drift vector 7 = —%hp (see text 
below eq. The relation to the Fokker-Planck equation can be further elucidated 

by noting that the "vector potential" A = £VV is a complete gradient and thus can 
be "gauged" away by means of the transformation W 
form of dHJ): 



e tv W\ yielding an alternative 



dW h 



2 r 



w - VW'. 



(16) 



at 2 

If we interpret W(r,p; t) as a one-time probability density and introduce the probability 
current (sometimes called the probability flux in the literature) ^3] 



h 2 

J' = -ihpW' - yVtr, 

we may write (|16|) in the form of a continuity equation: 
dW 



dt 



+ V • J' 



-VW' 



;i7) 



(18) 



where the potential term acts as a source term and violates the local conservation law 
associated with the Fokker-Planck equation. 

The recurrence relation (JT5|l clearly shows that the expansion in H introduced by 
Kirkwood ^2] is also an expansion in powers of the gradient operator as emphasized by 
Uhlenbeck and Beth ^3]. In the general case one needs to expand in powers of fi or, 
equivalently, in powers of the gradient operator using fl"2|) or (|15jl . as we will do here. 
In special cases, however, the compact form of the equations (fT4*|) or (jlfi)) may be the 
starting point of other effective approximation schemes. 
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4. Testing the method of separation beyond the TF approximation for the 
x 2 y 2 potential 

We now want to check whether the method of the separation x- and y- motions in the 
channels Ej works beyond the TF approximation. In the TF approximation, where 
the integrand in Z(t) is simply exp(— tH), the boundary Q appears in the argument of 
a logarithm. Higher WK corrections to the TF term lead to a power-like dependence 
on Q, and it is important to confirm that all dependence on Q is cancelled in the final 
answer, at least for the leading terms in the parameter tQ 4 ^> 1. Here we consider 
this problem up to the second-order corrections. In the later sections we will consider 
corrections up to the order h 8 . 

As we stressed before, we apply the WK method to the motion in the central region 
(|^|) \y\ — Q) m hs full scope, whereas for the motion in the channels (\x\, \y\ > Q) 
this method will be used only to study the quantum effects on the "slow" longitudinal 
motion, which is adiabatically separated from the quantized oscillatory motion in the 
transverse coordinate. 

Integrating (fl~2j) with respect to t, we have: 



W 1 = - -p-VV 

t 2 

W >=2 



(19) 



with V(x, y) = \g 2 x 2 y 2 . Integration over p x ,p y and making use of the symmetry of the 
Hamiltonian with respect to the interchange x <-> y, we find that the contribution from 
Wi vanishes and 



/( i 1 1 

dT Wo e 



-tv _ nt 9 



I + t92 I 

"-'10 T ~2 21 



(20) 



3 

where we introduced the abbreviations 

dT = dxdydp x dp y (21) 
and (for m > n)§: 

I mn = 4 [ Q dx f Q dyx 2m y 2n e- l ^l 2 . (22) 



io jo 

These integrals can be evaluated straightforwardly after the substitution x = w and 
y = J2/t(u/gw) and with the help of the condition for the validity of the adiabatic 
approximation Qt 1 ^ ^> 1. Note that this inequality does not contradict the condition 
permitting the use of the Wigner representation TiQt <C 1 if M 3 ^ 1. We obtain 



no » -^Q\ In - ^Q 2 . (23) 

§ Note that the case m — n needs to be calculated separately from the case m > n for the leading 
terms; see below. 
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and finally for Z 2 (t) in the square: 



-K-(gMQ) 2 . 



(24) 



For the sake of completeness and later use, we note the systematic structure of the 
integrals I mn . For a given value of m — n there are m — n + 1 such expressions: 
I m -n,o, I m -n+i,i, ■ ■ ■ , h(m-n),m-n- Applying the same method as above yields the 
following result for these expressions at fixed m — n: 

(2n - I) 



m 



n 



-Q 



2(m—n) 



(25) 



(gt l / 2 ) 2n + l 

For the second region (the four channels \x\ € [Q,oo] or \y\ G [Q, oo]), the 
integration is more involved. WK corrections to the TF term introduce pre-exponential 
functions of p x ,p y ,x,y and t in W 2 - Because of the fourfold symmetry of the Hamiltonian 
(J2J) it is sufficient to consider the channel x > Q. Closely following the method used in 
IB] we consider the motion in the x variable as free and the motion in the y variable 
as that of a harmonic oscillator with an x-dependent frequency with the Hamiltonian 



i i 

2 i 2 2 

^Py + ^xV 



(26) 



with oj x = gx and eigenvalues e n (x) = (n + \)hgx. 

In order to reduce the integrals over p 2 y and y 2 to derivatives of the well-known 
exponentiated sum-rule for the harmonic oscillator, we rescale the kinetic energy term 
in H y by a factor b and the potential energy term by a factor a. We denote the rescaled 
Hamiltonians by H' y . We then can generate any pre-exponential powers of p 2 and y 2 
by differentiating with respect to a and h and setting a = b = 1 in the final result. 
Integrating W% from (JT§)l over p x and x and performing the described substitutions and 
differentiations we obtain : 



dp x 



dx W 2 Tr (e 



-tH> 



dxD(x) 



(27) 



with 



D(x) 



x 

T 



2x 2 d 
3 da 



+ 



d 2 



3g 2 tx 2 da 2 
+ x 2 



d 2 



Tr (e 



a=b=l 



(28) 



3 db dadb 

In f)27jl we may discard the third term as it is 0(1/ xH) 1 with respect to other terms 
[x > Q, QH ^> 1 in the channel). || We note that the algebraic trick dealing with terms 
including the "fast" variables y, p y effectively reduce the problem in the channels to the 
level of TF terms with the rescaled Hamiltonians H' y . Now we can use the exponentiated 

|| One may notice that contributions arising from the derivative with respect to the "fast" variable y 
in Wi are not small, whereas derivatives with respect to x are kinematically negligible in the channel 
where Q 4 t ^S> 1. 
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sum-rule for the harmonic oscillator (see [7j) to obtain 

- 4 2a; 2 d 1 



D(x) 



x 



4 sinh £ 



3 da sinh (^^/a) 
x 2 d 2 1 



(29) 

2 dadbsmh^Va^)\ a=b=1 ' 
where £ = hgtx/2. Performing the differentiations, setting a = b = 1, and multiplying 
by 4 to account for all four channels, we obtain for the channel contribution to ^(t): 



Z 2 (t) = AK / edt, 



1 



11 cosh£ 



sinh £ 6 sinh 2 £ 
+ £ 2 £ 2 cosh£ 



-r— . (30) 
2sinh£ sinh 3 £j' V ; 

with £ = hgtQ/2(<^ 1). Performing the integrations in (|3Uj) (see (TB], integrals 2.477.3 
3.523.1, 2.479.4, 2.477.1, 3.523.1), we obtain for the channel contribution to Z2: 



Z 2 (t) = K 



-(^tg) 2 -2ic(3) 



(31) 

and (|3*T|) one sees 



where ^(z) is the Riemann zeta function (C(3) ~ 1.202). From 
that the Q-dependence of Z 2 (t) cancels, and the final answer is 

Z 2 (t) w -25.2K. (32) 

We achieved the Q-independence of the second-order correction to the partition function 
of x 2 y 2 model using the method of [Zl [Hj - The huge renormalization of the TF term in 
looks suspicious. 



5. Improved quantum motion in the channels 

The unexpectedly large coefficient in (|3*T|) suggests that we need to find a better 
treatment for the motion in the channel. The previous study of the quantum motion in 
the channel neglected the quantum nature of the "slow" motion along the channel axis, 
which was assumed to be free in [3 IE] • We already stressed in the Introduction that this 
motion is by no means free, rather, it is influenced by an effective linear potential created 
by the quantum fluctuations in the direction(s) orthogonal to the channel axis. Indeed, 
in the region \x\ 3> \y\, where the derivatives with respect to x are small relative to the 
ones with respect to y, we may first average the motion over the quantum fluctuations 
of y [E| described by the Hamiltonian H y (I2fij) and by the corresponding wave function 

My) = TPS (^f 1 e-^H n (y^Tn), (33) 
where H n (z) are the Hermite polynomials. The corresponding average value of H y 

(n\H y \n) = [n H — )%gx (34) 
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then becomes an effective potential for the description of the motion in the "slow" 
variable x: 

( h 2 d 2 1 \ 

+ ( n + n)h9x 4>n{x) = E(j) n {x). (35) 



V 2 dx 2 y 2 

This is the well-known Schrodinger equation for a linear potential with solutions in 
terms of Airy functions. Equation (|35|) clearly shows that, quantum mechanically, the 
"particle" is linearly confined along the channel axis and its motion is not free.^f The 
eigenvalue problem (J55J1 has a discrete spectrum. Note that this argument constitutes 
a sixth proof, in addition to the five proofs listed in [S], for the discreteness of the 
spectrum of the Hamiltonian (JHJ). It explains the large number in ()31)1 and (}32|) as an 
artefact of the assumption of free motion in the x-direction. This assumption becomes 
increasingly poor for the higher-order quantum corrections, leading to poor convergence 
or even divergence of the expansion in powers of h. Treating the "slow" x-motion in 
the channel adiabatically, we apply the WK expansion and eq. (112}) to the effective 
Hamiltonian 

Ht ] = \pl + {n+ l -)hgx (36) 

describing the motion of the nth quantum mode in the region x > Q. For the partition 
function of the one- dimensional motion of the nth mode we have (denoting p x simply 
by p): 

I P OO P oo 

Z( n \t) = — dp dxe-^ p2t - (n+1 2 )n9Xt W{x,p;t). (37) 



2nh 

Expanding W(x,p;t) in powers of % for the linear potential in (}36j) and using eq. (fl2]l. 
we easily obtain: 

W 2 = ^a 2 n t 3 (4-3p 2 t); (38) 

W 6 = 2 io. 3 1 3.5 !! a " t9(320 _ 720pH + 180pH2 ~ 9pH3) > 

W 8 = ^ 4 7n a^ 12 (8960 - 26880p 2 t + 10080p 4 t 2 - 1008p 6 t 3 + 27p 8 t 4 ), 

where a n = (n + \)Tig. The terms with odd indices contain odd powers of p and give 
vanishing contributions after integration over p. Carrying out the integration over p in 
(13 7j) . summing over n, and retaining all terms with derivatives up to order h 8 , we obtain 

Note that the phenomenon called "confinement" here is not the same as the phenomenon commonly 
referred to as quark confinement. In our case, the potential depends linearly on the field amplitude x(t), 
not on a spatial coordinate. We already emphasized the dual role of the coordinates x, y, z earlier, but 
we stress this point again here to avoid misunderstandings. One may also have called the phenomenon 
discussed here "self-confinement", as the fields themselves "prepare" the effective potential barrier 
prohibiting the escape to infinity. 
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for the contribution from a single channel: 

„ , N r°° dx ( h 2 t d 2 hH 2 d A , . 

ft 6 t 3 <9 6 ft 8 * 4 «9 8 \ 1 

+ n g + n 5 + 



82944 9a; 6 7962624 <9x 8 / 2sinh£' 

where again £ = hgxt/2. The expansion parameter h 2 t/Q 2 <C 1 in ()39)) is the squared 
ratio of the "diffusion" length mentioned in Section 3 and the separation scale Q. 
After performing the differentiations and integrating over x we finally have for all four 
channels: 



Z(t) = K 



t ( \2 \4 \6 

4incothi - __& + _ — a 3 + — — a 5 




(40) 



2 \2 3 -3~ ? ' 2 9 -3 2 ? 2 14 -3 4 ? 

1 

sinh^ 

with u = hgtQ/2 and A 2 = g 2 hH 3 . This expression, with the exception of the logarithm, 
can be written as a power series in u by substituting the Bernouilli expansion of the 
hyperbolic cosecans: 

-r -r oo r>2n— 1 I 

-T7 = 7 - 2 E Vt^^ 2 - 1 . (41) 
smh^ i (2n)! 

where i? 2n denotes the Bernoulli numbers (B 2 = 1/6, = —1/30,B 6 = 1/42, B 8 = 
—1/30,...). It is then easy to see that there are three kinds of terms in the series 
involving derivatives <9| n_1 in (J40j) : 

(i) terms with a positive power of u 2 ; there is an infinite number of such terms; 

(ii) terms with an inverse power of u 2 ; there is one such term at for each derivative; 

(iii) terms independent of u, and thus independent of Q; again there is only one such 
term for each derivative. 

The terms of type (i) are of the order \ 2n (hQt) k ~ 4n , where k denotes the overall power 
of h. The leading term in Z(t) at the same power of h being of order (hQt) k , arising 
from the expansion of the logarithmic term in (|40p. the contribution from the series 
is suppressed by a factor {tQ i )^ n 1. For the terms of type (ii) one can show that 
they are of the order (hQt) rn+1 (h 4 t 3 ) n . For a given power k = An — m — 1 of h their 
ratio to the leading contribution {TiQt) k is again of order (tQ A )~ n 1. Finally, for 
the terms of the kind (iii), their ratio to the leading terms at a given power k = An is 
(hH 3 ) n (hQt)- n = {tQ A )~ n = {tQ A )- k/A < 1. We thus conclude that each term arising 
from the series of derivatives in (|4()jl is smaller than the contribution from the leading 
term, and the relative suppression increases with k or n. However, it would be wrong to 
conclude that these subdominant terms can all be neglected, because there are certain 
terms among those of type (iii), which are independent of the cutoff Q and which are 
not canceled by similar contributions from the channels. 
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Having established that the terms arising from the logarithmic term in (|4U|) 
dominate in the expansion of Z(t) in powers of h, we now proceed to give those explicitly. 
The integrated form of the series (|4*T|) is: 

lncoth |.- ln | + gP^ B2 ^, (42, 

This series converges for u 2 < it 2 or (hgtQ) < 2tc, which is consistent with the condition 
permitting the use of the Wigner representation in the square [-Q, Q] and with the 
inequality tQ 4 ^> 1. We thus get for the partition function inside the four channels, up 
to terms involving powers of (tQ 4 ) -1 <C 1: 

4hi^-r + — (hgtQ) 2 - — — - — r:(hgtQ) 4 (43) 
hgtQ 12 v y * J 2 7 -3-5!P y W y ' 

31 -( h 9tQ) 6 - n J 2 J n J hgtQ) 8 + ■ ■ 



^0+2+4+6+8 (t) — K 



2 9 -3-9!P ^ ^ 2 16 -5-9!! 
We retained in ()43|) the higher-order terms for future use. 

The contributions Zq and Z 2 for the central square are given by (f2~4*j) and by eq. (5) 
from 7j. Their sum is: 



z{r + f Q \t) = K 



2 In V2t l ' 2 g 2 Q 2 + C - ^{KgtQ? 



(44) 



Adding the first two terms in P3*|) and (j44j) we finally obtain for Z up to the order of h 2 
the expression (J3J) found in j7j, i. e. the Q-dependent corrections to the TF term vanish 
at the second-order of h. Below we shall see that this statement is correct even up to 
the order h 8 . We conclude that the unusually large number renormalizing the TF term 
in (J32J) is an artefact of the neglect of the quantum character of the motion along the 
axis in the hyperbolic channels. 

6. The corrections to the TF term vanish up to order h 8 

In this Section we show that the leading Q-dependent quantum corrections to the TF 
term up to the order h 8 cancel if we treat the quantum mechanical motion in the channels 
correctly. Actually, we already gave the corrections to Z(t) up to h 8 for the correct 
motion in the hyperbola channels (see formula ()43j) ) . It remains for us to calculate 
these corrections for the square x,y G [— Q, Q], a straightforward though cumbersome 
task. First of all,we need to know W±(f) using ()12|) . There are two types of contribution 
to W&{t): With m — n = 2 and m = n. Terms with m — n = 2 give the main contribution 
of the order tQ 4 , terms with m = n contain only logarithms of tQ 4 and may be neglected 
with the precision \ia(tQ 4 )/tQ 4 1. Here we give the expression for W^it) after the 
integration over p x and p y : 

2ng 2 t 3 fl 2 11 



l drWie -v, = ^±(l 9hIw .^ igH)%l 



+ ^(g 2 t) 3 h 2 -^l 00 + g 2 tl u 

-g(^) 2 / 22 + ^(^) 3 /3 3 ) (45) 
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Integration of /„ 
we obtain: 



'oo 



over x and y for m 7^ n was already done before (see ([25])): for m = n 

/ 27T r 

" ln(/Q 4 t) + C + ln2 



hi 



W 2 



ln(/Q 4 t) + C + ln2-2 



/OO 



'33 



0^)5/2 

WF 2 



ln(g 2 QH) + C + ln2 



8 



46 



ln(^g 4 t) + C + ln2- — 

15 



(46) 



Collecting all terms together we have from (|45|) for the corrections to Z(t) of the order 
of h in the square x,y E [-Q, +Q]'- 



},Q] 



(t) = K 



7 



g^tQ" 



(47) 



_2 7 -3-5!!- 

+g 2 hH z [\n(g 2 Q 2 -8(7- 16 In 2] 

where the first term comes from I mn with m — n = 2, and the logarithmic terms, as 
we remarked before, arise from the contributions I mm . We see that the second term in 
([47]) is of order In^tQ 4 ) / Q 4 t <C 1. Discarding it and adding the h 4 correction from the 
channels in ([43)1 we find that with the precision ln(tQ 4 )/tQ 4 <C 1, 



-Q,+Q] 



+ Z 



[Q,oo] 







(48) 



as it was for the second-order corrections (with the precision l/tQ A ). 

Consider now the /i 6 -order corrections to the partition function. In ([43]) we included 
the dominant contribution from the channels [Q, 00] (fourth term). For We in the square 
we have, after integration over p x and p y and using the notation of 

2vr(^) 6 r „ , 249 



dT Wa e 



-tv 



t-9\\ -2 6 
119 



-6 1/30 + 



<9 2 t)l 



11 



(^) 2 / 5 2 + ^(^) 3 / 6 3 



(49) 



35 

4 w 7 " u " ' 24 ( 

Using ([23]) for the evalulation of the I mn and collecting all terms, we obtain for Zq~ q ' q ^ 
in the square: 



-Q,Q] 



(t) = -K- 



31 



-{hgtQf. 



J 6 V) — _J1 3 . 2 9 . 9!! vv^J ■ (50) 
Once more, the contribution from the square is exactly canceled against the one from 
the four channels given in ([43)1 . as it happened for the second and fourth order: 



zlr Q ' Q] (t) 



zf>°°\t) 







(51) 



with the precision ln(tQ 4 )/(tQ 4 ) < 1. 

Finally, we give here the order h 8 corrections to z\ in the square: 



dT W* e 



-tv 



2n(gtf 
t-5\\-2 7 



1296 



259 2 893 
140 _ 540^ thl + 5040 



1261 
7560 

23 (9 2 t) 3 hs+ ° 



(g 2 t) 2 h2 



10365 



(52) 
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where we again discarded terms proportional to I mm , which are of the order 
ln(tQ 4 ) / (tQ 4 ) <C 1 with respect to the terms containing I mn with m > n. Using now 
the general expression for I mn we obtain: 

/ Jy "-"° iiroi W' (53) 

and get the result 

1 97 

Z [ f Q ' Q] (t) = K h 21% m (hgtQf (54) 

which miraculously cancels with the contribution from the channels in ([43)1 . 

Although we cannot prove such a cancellation in general, to all orders, we have no 
doubt that all higher-order corrections to the partition function containing powers of 
(hgtQ) cancel in the limit tQ 4 ^> 1. Of course, there are other corrections involving 
powers of (tQ 4 )' 1 , but these are suppressed due to the classical condition of adiabaticity 
tQ 4 3> 1. We shall consider this issue in Section 7. 

In anticipation of later applications, we remark here that there is a strong 
correlation between the power of H, denoted by k > 2, and power of the dominant 
terms in the limit tQ 4 3> 1 . A systematic analysis of the higher-order corrections using 
Mathematica leads to the conclusion that the difference between m and n in the leading 
integral I mn from (J22|) is m — n = \k. Besides the terms with the largest difference 
m — n, which give the leading contribution to Wk and which we retain in (J49|) and (J52)) . 
there are also terms involving I mn with m — n = |& — 2£ with £ = 1, 2, ... < ~k. There 
is also a correlation between the powers of t and m, n. The analysis shows that for the 
terms with m — n = the power of t is 2m — n — 1, and for terms with m — n = ^k — 2£ 
the power is 2m — n — 1 — 3£. The factor of g for I mn is g 2m . Straightforward calculations 
similar to the above show that the ratio of the contributions to the partition function 
for the terms m — n = |A; — 2i to the ones for m — n = k/2 is of the order 1/ (tQ 4 ) e 1. 
Indeed, we find: 

Zt' n \t) = KihgtQf^^ (m-n = ^- (55) 

Zt' n \t) = K(g 4 tQ 4 y\hgtQ) k ^^ (m - n = \ - 2*) . 

This justifies neglecting these subdominant terms in our calculations here. 

Finally, at any power k of h terms with logarithms (m = n) may be discarded with 
respect to the terms with powers of Q 2 . Indeed, the last terms with m > n are of the 
order of (hgtQ) k , terms with m = n are of the order of (g 2 hH 3 ) k/4 ln(tg 4 Q 4 ) and may 
be neglected with the precision \ia(g 4 tQ 4 ) / (g 2 tQ 4 ) 1. Again, with increasing k the 
precision improves. 

7. Q-independent terms in the channels and central region 

We already mentioned that the derivative expansion for the channel contribution ()40|) 
has one Q- independent term for each power of of A 2 = g 2 h 4 t 3 . From the expansion (|41j) 
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it is evident that these terms have the form 

2 2n-l _ l 



n(2n V . P"- 1 )"^"- <«> 

Each terms is of the order (tQ 4 )^ 1 <C 1 compared with the leading term from the 
expansion of the logarithic term in (140(1 . Collecting these terms we obtain a series of 
Q-independent terms contributing to Z(t): 

1 2 3 — 1 2 5 — 1 2 7 — 1 

- — B 2 X 2 + — — -5 4 A 4 + -u — --BgA 6 + — — -B 8 \ 8 + ■■■ (57) 
2 3 • 3 2 -3 2 • 3 5 2 • 3 5 

We can rewrite this series in the form 

^ 2 2 "(2 2 - 1 -l)(2n-l)!! D / A \ 2n 

h 2^)n{2n)\ ^{Ws) ' (58) 

The factor (2n — 1)!!, which arises from the (2n — l)-th derivative in (|4"Uj) . spoils the 
convergence of this series for Z(t). Using the identity 2 n n\(2n — 1)!! = (2n)\, we obtain 
for the series of Q-independent terms: 

Z [QM = K f^ 4(2 2 - 1 -l) / A \ 2n 

This series is only asymptotic despite the smallness of A 2 , because the Bernouilli numbers 
grow factorially at large n: 

2{2n)\ 

W) 

Next we turn to the Q-independent contributions to Z(t) from the central 
region. We start with the investigation of the corrections to the dominant terms. A 
straightforward calculation yields the following correction factor to I mn from (J25j) : 

1 " w^W is ' QH) " (61) 

This factor plays an important part making the corrected I mn well behaved in the limit 
m = n. Indeed, we note that with e = 2{m — n): 

(2m -1)!! ™ 1 

the diagonal elements (m = n) become 

/ 2vr(2m - 1) 



b» - Ssr- (60) 



lim I mn -— f 

e-> o mn (gH) m+ ^ 



2 



h 9 2 qh)-2y: 



i 2^-1. 



(63) 



demonstrating that the corrections to I mm are independent of Q, with the exception of 
the logarithmic term. One easily confirms this general result by explicit calculations. 
Substituting x = w and y = J2/t(u/gw) in (f2"2"j) we get: 



2m + l 

2 \ 2 rQ dw r w 



Imm = 4 ( ] / — / duu 2m e- u \ (64) 
g HI Jo w Jo 
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We may write this expression as a derivative of the error function: 



2m + l 

2 



w ° dw 



4(— ) ' ^(-1) 
\g 2 t 2 v 1 da m Jo ^/aw 



erf (y/aw] 



(65) 



0=1 



where Wq = g 2 QH/2. Integrating by parts, discarding the exponentially small terms in 
Wq, then differentiating with respect to a and finally setting a = 1, we obtain: 



2tt 



(g 2 ty 



2 m r°° 

(2m-l)!!lnw 2 - —= \ dtt m -^e' l \nt 

■J TV JO 



(66) 



For the integral in the second term of this expression one finds (see ref. [TH] formula 
4.352.2): 

1 



dt t m ~^e~ l hit 



(2m - 1)!! 



2E 



2£-l 



C-21n2 



which yields 



2vr(2m- 1)!! 



ln(/Q 4 t) +C + ln2-2^ 



21 



(67) 



(6J 



Next we use a trick, first introduced by Euler, replacing the finite sum by an asymptotic 
series (see e.g. [T9]): 



E 



21 



1 

T ~ 2 



oo q2€-1 



C + ln(2m) + E 



- 1 



5m 



2V 



2/ 



(69) 



We are now ready to combine the contribution from the central region with the 
asymptotic series obtained earlier for the Q-independent contribution for the 
channels. A special case is the case m = n = 0, for which there is no infinite sum, 



giving 



i - r 



\n(g 2 QH) + C + \n2 



(70) 



Taken together with the logarithmic term from the channels, this leads to the Q- 
independent expression (J2J) obtained in the improved TF approximation Our detailed 
analysis has shown that the structures I mm containing Q-independent terms appear only 
in with k — 1,2,.... This means that such terms appear only at the orders h Ak and 
implies that the expansion parameter of Z(t) is A 2 = g 2 h t 3 . 

For a given power A 2 ™ there are 3n quantities I mm (m — 1,2, ... , 3n) and I 00 , which 
we collect separately. The results given above allow us to write the Q-independent 
asymptotic series from the central region. Including also the TF term (j3J), we thus 
obtain: 

1 



Z [ - Q ' Q] (t) = K 



In 



+ 91n2 + C 



g 2 h 4 t 3 

oo 

+ E x2n (4 n) (C + ln2) 



n=l 
3n 



m=l 



oo o2£-l 



lnm + 2J 



6m 



- 1 

2V"" 



-B 



2C 



(71) 
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Here the a$ (m = 0,1,2,..., 3n) denote the numerical coefficients in (g 2 t) m I mm 
occurring in Z 4n . For example, for n = 1 we have four such numbers with alternating 
signs: 

a ° ~ 60' ai "16' ° 2 ~ 720' fl3 " 576- 
Summing up all these results, we may surmise that in the limit tQ 4 ^> 1 all Q- 
dependence is canceled, and we are left with two asymptotic series (fSUjl and (J7T|) from 
the two regions contributing together to the partition function of the two-dimensional 
YMQM model. We have shown explicitly how the small parameter (tQ 4 )^ 1 , artificially 
introduced in the adiabatic separation of the degrees of freedom in the channels, 
transmutes into the genuine quantum mechanical parameter A 2 = g 2 h 4 t 3 . 



8. Conclusions 

We have shown that the richness of the classical YM mechanics with a x 2 y 2 potential 
translates into, and even gets amplified by, the quantum mechanical properties of the 
system. The YM quantum mechanics exhibits a confinement property, which strongly 
influences the quantum mechanical motion in the x 2 y 2 potential. At higher order in % (up 
to h 8 ) this results in the vanishing of the leading quantum corrections (for tQ 4 1), 
when we correctly take into account this property for the motion in the hyperbolic 
channels. We believe that this result may survive to even higher order. We also derived 
a novel form of the equation for the Uhlenbeck-Beth function W(r,p;t), which is the 
basis of the WK expansion. We believe that our expression can be a starting point 
for new approximation schemes using techniques from diffusion theory. We hope that 
the lessons derived from the present study of the higher-order quantum corrections 
to the homogeneous limit of the Yang-Mills equations will be useful for an improved 
understanding of the internal dynamics of the Yang-Mills quantum field theory. 
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